Hall resistance and Lorentz ratio models in YBa2Cu307 
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Two-dimensional models of heat capacity, conductivity («:), Hall resistance {Rh) and the Lorentz 
ratio (£) have been derived using the ionization energy {Ei) based Fermi-Dirac statistics (iFDS) 
for overdoped Cuprates. These models reproduce the experimentally measured p{T) and Rh{T). 
The variation of C is weakly T-dependent due to the experimental k.(T). The e-ph coupling in 
the electrical resistivity has the polaronic effect that complies with iFDS, rather than the e-ph 
scattering, which satisfies the Bloch-Griineisen formula. These models are purely Fermi liquid and 
are not associated with any microscopic theories of high- Tc superconductors. 
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1. Introduction 

The electrical properties of higli-Tc Cuprate supercon- 
ductors (HTSC) are intrinsically enigmatic in both ex- 
perimental and theoretical framework compared to other 
oxide compounds, including manganites. Partly due to 
its huge potential in a wide variety of applications, in- 
tense focus is given on the nature of conductivity of these 
materials to shed some light on the puzzling temperature- 
dependence issues in heat capacity (C), heat conductiv- 
ity, resistivity, Hall resistance and Lorentz ratio. The 
conflict in term of T-dependency between p and Rh 
is also one of the unsolved mystery in HTSC. Even 
though this paper does not solve it microscopically, but 
it points out that the Fermi liquid with strong correla- 
tions is not downright incorrect, at least for over- and 
optimally doped HTSC. Apart from HTSC, the applica- 
tions of ionization energy {Ej) based Fermi-Dirac statis- 
tics (iFDS) on ferromagnets, diluted ferromagnetic semi- 
conductors and doped ferroelectrics have been derived 
and discussed analytically II II mill II QL The iFDS 
in HTSC have been successfully exploited [ll, M jdl via 
the experimental data reported in the Refs. la, S 0, 
miHIIllIllIllISIIllIllIi- Recently, the said 
puzzling T dependencies as well as the spin gap phe- 
nomenon have been tackled with the coupling of iFDS 

and charge spin separation mil 113. Unfortunately, 

the pure charge-spin separation, is believed to have se- 
rious shortcomings [2I HI HI H^. In addition, 
there are also numerous experimental reports with con- 
troversial interpretations surrounding the spin Pseudo- 
gap phenomenon 0,S0|3Ol|3ll|3^|3i|34- As 
such, by ignoring the spin Pseudogap phase, the thermo- 
magneto-electronic properties of YBa2Cu307 will be dis- 
cussed based on the iFDS by heavily relying on the basic 
transport experiments such as the resistivity, Hall resis- 
tance, heat capacity and heat conductivity. It is inter- 
esting to note that these purely Fermi-liquid models are 
able to reproduce the related experimental data reason- 



ably well even if they are only for over- and optimally 
doped HTSC. The polaronic effect that arises as a re- 
sult of iFDS is solely due to heavier effective mass effect, 
which could indicate the existence of polarons. But this 
indication is just an extrapolated assumption since heavy 
electrons do not necessarily form polarons. 



2. Theoretical details 

The free-particle Hamiltonian of mass m moving in 3- 
dimensions is given by 
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Here, we have make use of the linear momentum oper- 
ator, p = — iftV. Subsequently, one can write the time- 
independent Schrodinger equation for the same particle, 
however in an unknown potential, V{x) as 
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In the second line of Eq. (O, one can notice that the 
influence of the potential energy on the total energy of 
that particular particle has been conveniently parameter- 
ized as ^. This energy function, ^ will be characterized 
later in such a way that one can replace E + V{x) with 
Eo+ (, in which, Eq ^ E at T ^ 0. Add to that, from 
Eq. ||2Jl, it is obvious that the magnitude of ^ is given by 
±^ — E ~ Eq + V{x). Physically, it implies the energy 
needed to overcome the potential energy as well as the 
bound state. Literally, this is exactly what we need to 
know in any condensed matter, i.e., this magnitude is the 
one that actually or reasonably defines the properties of 
the quasiparticles. Subsequently, we obtain 
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fc^ = {2m/h'^)[Eo ± i? and i5o in a given system 
range from +oo to for electrons and to — oo for holes 
that eventually explains the ± sign in ^. Now, Eq. Q 
can be solved to give 

if = C'n X 

exp[i{ko,^ ± k^.a;)x + i{ko,y ± k^,y)y + i(A;o,2 ± k^,z)z]- 
^k(o.,) = C.e^^'^^y^. (5) 

fc2 = (kl, ± klJ + (kly ± kly) + (fco% ± kl,). By 
employing the orthonormality and Plancherel's theorem, 
one can find the normalization constant, Cat by compar- 
ing Eqs. © and as shown below. 

(■^ko I '^ko±k,) = / / / v'ko * '^koik.'^^^y^^ 
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e'^° '^^-^">dk^dkydk^ = S{r - r'). (7) 



Hence, Cn — l/(27r)'^/^. finally, the normalized wave 
function, which corresponds to Eq. (O is 
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In a physical sense as stated earlier, = E — Eo + V{x), 
is in an identical scale with the energy needed to free an 
electron from an atom in a given crystal. As such, we 
apply the concept of ionization energy where, ^ = jrjreai 
— Ej + V{x), to justify that an electron to occupy a 
higher state N from initial state M is more probable 
than from initial state L if condition Ei{M) < Ej{L) at 
certain T is satisfied. As for a hole to occupy a lower 
state M from initial state N is more probable than to 
occupy state L if the same condition above is satisfied. 
It is well known that the exact values of Ej are known 
for an isolated atom. In this case (for an isolated atom). 
El can be evaluated with 



However, substituting the same atom in a crystal gives 
rise to the influence of V{x) and in reality, iJJ'^'*' cannot 
be evaluated from Eq. 0. Nevertheless, the of an 

atom or ion in a crystal is proportional to the isolated 
atom and/or ion's Ei as written below. 
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It is this property that enables one to predict the vari- 
ation of electronic properties of superconductors with 
substitution reasonably well. The constant of propor- 
tionality, a is a function of averaged V{x) and varies 
with different background atoms. For example, in 
YBai_a;Caj,Cu2 07 system, YBai_a;-Cu207 defines the 
background atoms or ions. Therefore, one needs to em- 
ploy the experimental data to determine the magnitude 
of 5 = iJJ'^''' = Ei + V{x). 

Recall that Eq. Q simply implies that the one-particle 
energies Ei, E2, Em for the corresponding one- 
particle quantum states qi, (72, Qm can be rewritten 
as {Eo ± Ei)i, {Eo ± Ei)2, (£^0 ± £^/)m. It is also 
important to note that Eq + Ej = E^iectrons and Eq — 
Ej = Ehoies- As such, for n particles, the total num- 
ber of particles and its energies are conserved and the 
conditions to fulfill those conservations are given by 
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dui = 0. 
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J2iEo±Ei),n,^E, Y.(^o±Ei),dn,^0. 

i i 

(12) 

Subsequently, the Fermi-Dirac statistics based on ion- 
ization energy can be derived as 
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(13) 



By utilizing Eq. (jTSl and taking cxp[^ + A(i?±£;/)] > 
1, one can arrive at the probability function for electrons 
in an explicit form as 



/e exp 



— yU — A 



21,2 



2™ 



E, 



(14) 



Similarly, the probability function for the holes is given 



by 
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fh = exp 
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The parameters n and A are the Lagrange multiphers. 
h = h/2-K, h = Planck constant and m is the charge 
carriers' mass. In the standard FDS, Eqs. (|14|1 and l|15|) 
are simply given by, /e(k) = exp[— — \{h?'k'^ /2m)] and 
//j(k) = exp[^ + A(/i^k^/2m)] respectively. 

Subsequently Eq. can be rewritten by employ- 

ing the 2D density of states' (DOS) derivative, dn — 
i^ko(iko/27r, Eqs. ((njl and ifTsl) . that eventually give 
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fe{E,Ei) = exp 



A (£,£;,) = exp 



Ep — Ej — E 



E-E, 
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Note that Eqs. (^0)) and if^ simply imply that 
lie{iFDS) = Ai(r = 0) + \Ei and fih{iFDS) = n{T = 0) 
— Ai?/. In fact, fi{FDS) need to be varied accordingly 
with doping, on the other hand, iFDS captures the same 
variation due to doping with XEi in which, fi{T = 0) is 
fixed to be a constant (independent of T and doping). 
Furthermore, using Eqs. H12|l and H22|l. one can ob- 
tain 



Note here that Eq is substituted with /i^ko/2r7i. n and ^ 
p are the respective concentrations of electrons and holes. 

denotes area in k-space. The respective solutions of 
Eqs. (fTB)l and lfT7|) are given below 
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Equations (fH^ and itT^ respectively imply that 
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Eventually, Eq. H25|l proves that the total energy of n 
particles considered in both FDS and iFDS is exactly the 
same. 



(20) 



3. Heat capacity and its conductivity 



^lH{iFDS)^ ^lh~XEI. (21) 
On the other hand, using Eq. p2(l . one can obtain 
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Equation 122|l . after appropriate substitution with 
Eq. p8l) is compared with the energy of a 2D ideal gas, 
E = nksT. Quantitative comparison will eventually 
leads to XiPDS = Xfds = ^/ksT where /cb is the Boltz- 
mann constant. The distribution function for electrons 
and holes can be written explicitly by first denoting fi — 
—Ep (Fermi level), A = l/ksT and substituting these 
into Eqs. (|14|l and H15|) will lead one to write 



Electrons and phonons can be excited to a higher en- 
ergy levels satisfying the ionization energy based Fermi- 
Dirac {fiFDs{E)) and Bose-Einstein fsEsiE) statistics 
respectively. Therefore, the heat capacity can be explic- 
itly written as 
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(26) 



h = h/2'K, h denotes Planck constant, while m* rep- 
resents the effective mass. The respective distribution 
functions for EES and iFDS (using Eq. Q) are given by 
fBEs{E) = l/[exp [{E - EF)/kBT] - 1] « exp [{Ep - 
E)/kBT] emdf.FDsiE) - l/[exp [{E-EF+Ei)/kBT] + 
1] « exp [{Ep — El — EJ/ksT]. These approximations 
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are necessary to avoid the exponential integral function, 
'E.i{z) = J_^[e~* /t]dt, which has a branch cut discontinu- 
ity in the complex z plane running from ~oo to 0. Ad- 
ditionally, I highlight that for classical particles satisfy- 
ing the Maxwell-Boltzmann statistics (MBS), there is no 
such thing as Ej. Consequently, one should not assume 
that exp[^ + \{Einitiai state ± Ej)] > 1 should give the 
MB distribution function as a classical or a free-electron 
limit. One can indeed arrive at MBS by first considering 
the additional constraint, £^7 = in such cases, where 
Etotai now equals to E identical with the standard FDS 
and MBS. Therefore, the electron's Fermi level, Ep term 
that contained in fsEsiE) corresponds to the phonons' 
energies above this Ep, in which this Ep does not imply 
phonons' chemical potential. In other words, phonons 
with energies < Ep are neglected. In summary, iPDS 
captures the Fermi liquid 7^ 0) rather than the Fermi 
gas {V{x) = 0). 

The total heat capacity in Eq. (|26|l has been written 
as C = C*^ -I- C^'' as a result of the total heat current, 
« = -EaJQ(VT)-i = Ec.C"v^Fru/2. a = electron (e), 
phonon (ph) and — e-e, e-ph scattering, vp denotes 
the Fermi velocity and is the Boltzmann constant. 
Importantly, Ej is microscopically defined as 

eiO,k) = l+'^exp[X{E"p-Ej)]. (27) 

e(0, k) is the static dielectric function, k and A are the 
wavevector and Lagrange multiplier respectively. The 
Ep denotes the Fermi level at K, while the ICg rep- 
resents the Thomas-Fermi screening parameter. Unlike 
electrical resistivity in YBa2Cu307, its 2D heat conduc- 
tivity is equally strongly influenced by e-e and e-ph in- 
teractions, hence (after taking Ep = ^m*Vp) 



needs to be determined from other techniques, be it the- 
oretical or experimental. Therefore, only the electron's 
efl'ective mass is highlighted here. In other words, instead 
of addressing as the reduced mass of ions, it has 

been labelled as phonons effective mass so that one can 
conveniently identify it as the parameter belonging to the 
phonons' contribution. 

4. Resistivity and Hall resistance 

The equations of motion (EOM) for charge carriers in 
a6-planes under the influence of static magnetic (H) and 
electric fields (E) can be written in an identical fashion as 
given in Ref. |23 > which are given by m* [d/dt+1/ Te] ft = 
eEf, + eHc^Q and m*[d/dt + 1/Te\va — eE^ — eHcVh. 
The charge, e is defined as negative in the EOM above. 
Moreover, it is important to realize that the existence of 
electrons in a6-planes below Tcrossover are actually holes. 
The existence of holes in a6-planes was discussed inten- 
sively in the Refs. [13, 113, Hal- The subscripts a, b and 
c represent the axes in a, b and c directions while the 
subscript ab represents the a6-planes. In a steady state 
of a static H and E, dva/dt — dvb/dt — and Va — 
hence one can obtain E^ = eHcEbTe/m*. The Hall 
resistance and current along a- and 6-axes are respec- 
tively defined as i?^' = Ea/jtHc, jb ~ ^b/p in which, 
tan 6*^' = Ea/Efe. Parallel to this, i?^^ = tan6I^V/Hc. 
jb is the current due to holes motion along 6-axis and 
0^'' is the Hall angle in a6-planes. Furthermore, one can 
rewrite tan^^-* as tan6'^'' — ellc/m*AT'^, which eventu- 
ally suggests, cotS^-* (xT^. A is Te dependent constant 
and is independent of T. The 2D resistivity model, p{T) 
is given by (li, (2, l^j 
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(28) 



The explicit form of Eq. (|28ll can be obtained after 
substituting Eq. into Eq. (PSl appropriately. The 
electron-electron scattering rate, t~\ — t^^ = AT^ while 
the electron-phonon scattering rate, T~2pf^ is assumed to 
be proportional to T" in which a > 2. The Fermi- level 
in Eqs. (|26|l and H28|) implies that the phonons consid- 
ered here have the thermal energies in the order of or 
higher than the electrons' Fermi energy which eventually 
means that these electrons cannot form Fermi gas. In 
simple words, if the thermal energies of the phonons are 
less than the electrons' Ep, then these electrons can act 
as Fermi gas and one may employ the Debye approxi- 
mation. This is another reason why Debye model works 
extremely well at intermediate and low temperatures in 
common metals. However, the phonons' effective mass is 
equal to the ions reduced mass due to phonons interac- 
tion with free-electrons, i-/m*f^ = l/rTig + l/m.io„, which 



p(T) ^ A^Texp 



Ei + Ep 
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(29) 



Utilizing Eq. H29|) . one can show that the Hall resis- 
tance is given by 



Rh — 



m*TkBE 



exp 



Ei + Ep 



T 



(30) 



Thus, it is clear that Rh is proportional to 1/T re- 
gardless of the axes. Detailed analysis and diagnosis of 
Eq. (|29|l with a wide variety of experimental data are 
well documented in the Refs. Optimally doped 

YBa2Cu307_5 single crystal (Al) obtained from Ref. [33 
will be utilized in the following analysis. Equation l|29() 
has been employed to theoretically reproduce (indicated 
with a solid line in the inset of Fig. the experimen- 
tal Pab{T) by varying the T- independent scattering rate 
constant, A (7.3 x 10~^ ft-cm) whereas Ej + Ep = 
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Tcrossover (Tcr) IS taken as K, since any optimally doped 
YBa2Cu307 gives Tcr Tc in which Tcr is not observ- 
able from the resistivity measurements. I.e., Tcr cannot 
be predicted accurately from the normal state resistivity 
measurements. On the other hand, the ii^''-* (T) data and 
the plot using Eq. (|30|l are depicted in Fig. Note that 
A = A-Kh^/kBC^ from Eq. and Ah = TTh^/m*kBe 
= 347 JKCs^kg-i from Eq. In the latter approxi- 

mation, m* = 50mo, ttiq is the rest mass of the electron. 
In order to accurately fit the experimental R^^^'^ (T) data, 
the effective mass should be equal to 73mo, which in turn 
gives the charge carriers density as p = 8 x 10^^ cm^'^, 
in accordance with the Refs. |3l Il7|. 



5. Lorentz ratio 

Lanzara et al. have shown that the e-ph coupling 
is somewhat inevitable, which has been observed via 
ARPES technique. Indeed this supports the notion of 
polaronic effect above Tc in cuprates. One should note 
that the observation e-ph coupling does not mean that 
there is a e-ph scattering since normal state p{T) mea- 
surements thus far failed to reveal any e-ph scattering 
(strong T-dependence) . Actually, this is not because of 
p{T)'s blindness, but due to polaronic effect represented 
by Ej , which gives rise to the effective mass (m* ) of elec- 
trons instead of strong T-dependence. The heavier m* 
implies the existence of polaronic effect in the normal 
state of HTSC that also suppresses e-ph scattering but 
not the e-ph coupling in term of polaronic effect. Sim- 
ilarly, isotope effect (^^O, ^^^O) in cuprates [iollillli^ 
also reinforces the polaronic contribution via e-ph cou- 
pling rather than e-ph scattering. In fact, Hofer et al. |4(tI | 
claimed that m* reduces towards the optimally doped 
HTSC. This scenario is consistent with iFDS based mod- 
els that predicts Tcr also reduces towards optimal doping. 
Simply put, reduced Ej will eventually lead to reduced 
m* and consequently the influence of isotope doping on 
m* is less effective in optimally doped regime as com- 
pared to under doped. The inappropriateness of the e-ph 
scattering in YBCO7 will be discussed in detail based on 
the Bloch-Griineisen formula shortly. From the defini- 
tion, C can be written as 



C = —K = kA 



h^ 



T 



■ exp 



Ej + Ep 



T 



(31) 



Interestingly, Sutherland et al. [43j have reported only 
a slight increase (upward deviation) in afc-plane's heat 
conductivity with phonon contribution [n]^^ / 
< 1.3) above critical temperature (Tc) for overdoped 
YBCO. Their results will be used to discuss the accu- 
racy of Eqs. (|28|l and H31|) to capture the experimental 
data. 



6. e-ph scattering in resistivity 

Firstly, the Bloch-Griineisen (BG) formula will be re- 
visited in order to rule out the e-ph scattering in the 
normal state of YBCO7. Recall that the polaronic ef- 
fect that arises from the Ej based Fermi-Dirac statistics 
(iFDS) has been successful to explain and predict the 
evolution of resistivity with doping and to enumerate the 
minimum valence state of multivalent dopants in HTSC, 
ferromagnets and recently in doped- ferroelectrics. But 
iFDS does not reveal the inadequacy of the free e-ph 
scattering directly (only indirectly). Basically, accord- 
ing to the e-ph scattering, the electrons from Ba^"*" and 
Sr^+ as in Y(Bai_2,Sra;)2Cu307 has the same effect on 
transport measurements while iFDS points out that the 
kinetic energy (KE) of the electrons from Ba^+ is not 
equal with the KE of the electrons from Sr^+ , which gives 
rise to significant changes of resistivity with small dop- 
ing. Again, if one assumes KE (Ba^+) = KE (Sr^+), then 
the theory of the e-ph scattering is indeed applicable due 
to isotropy in KE (all the free electrons have an identical 
KE, which eventually defines the Fermi surfaces). Hence, 
to further evaluate the incompatibility of the e-ph scat- 
tering in YBCO7, the BG formula stated in Eq.|23 
is employed to plot the T-dependence of p{T) (assuming 
Tc^phi^D) cx Tc-ph{2D)) and CiT). 



Pbg 



1287rm*fcsT5 



eD/2T 



sinh X 



dx. (32) 



Xtr = electron-phonon coupling constant, m* — av- 
erage effective mass of the occupied carrier states, 6_d 
= Debye temperature, n = free electrons concentration. 
The C{T) can be simply written as 



C-BG 



KXt 



l2^iTm*kBT^ 



sinh^ 



-dx. (33) 



7. Analysis 

Figure [21 a) and b) depict the T-dependence of p{T) 
(Eq. ^) and C{T) respectively. The £(T) based on 
BG's approach after incorporating the experimental k are 
indicated with A (Od = 200 K), o (Od = 300 K), and □ 
{Qd = 350 K). On the other hand, the experimental and 
iFDS based theoretical plots (Eqs. if^ and ^ ) 

are shown with • and a solid line, respectively in Fig. 13 
b). Note that in Eq. H28|l . a = 3 is used complying with 
the earlier assumption of T"^^. This value is reasonable 
since a in the free e-ph scattering of conventional met- 
als are known to vary between 3 and 5, depending on 
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T's range that can be verified from Eq. H32|l . The ex- 
perimental C{T) is obtained from the resistivity [i^ and 
heat conductivity ^3 measurements of optimally doped 
YBCO. 

The inverse proportionality of the theoretically deter- 
mined K with T from Eq. H28|l is understandable since the 
electrical conductivity is proportional to 1/T and there 
are phonon contribution as well. As a result of this, C{T) 
is also inversely proportional to T. It is not possible 
to evaluate Eq. H28|l quantitatively due to the unknown 
magnitudes of Ej, Ep, m*i^ and A^-ph- However, the 
measured k in the normal state of YBCO hardly shows 
strong T dependence indicating the existence of some 
not-yet-known physical phenomena, which complicates 
our understanding of HTSC generally. Anyhow, by us- 
ing the experimentally determined k, one can verify the 
accuracy of the resistivity equations (between Eq. H29|l 
and (|32|l 'l. The former equation is entirely based on 
e-e scattering while the latter contains the essential e- 
ph scattering mechanism. To this end, the lorentz ratio 
based on iFDS (Eq. ^) and EG (Eq. ^) are com- 
puted using the almost T independent or experimental 

K. 

The iFDS model reproduces the T dependence trend, 
remarkably identical with the experimental data as op- 
posed to the BG's approach. Both iFDS and BG mod- 
els with the experimental k have been plotted in Fig. [3 
b), in which the latter model is plotted at different Qd- 
Eventually, one can convincingly state that e-ph scatter- 
ing mechanism is significantly negligible in the electrical 
resistance measurements. The plot that corresponds to 
Eq. with experimental k is obtained using Ej~\-Ep — 
10 K, (which is less than Tc as a result of optimal or over 
doping) and experimental n that eventually give A = \ x 
10~^ ri-cm. This magnitude is remarkably identical with 
the o ptim ally doped crystalline YBCO sample of Hagen 
et al. |l9| and Leridon et al. ji^ that have been calculated 
{AHagen,LerriAn-n = (1.1,1.4) X 10"* Vt-cvci) and reported 
in the Refs. 0,|21j. Importantly, even though Eq. ll^ 
can be shown to capture the experimental T-linear prop- 
erty of p{T), but it also fails to explain the Tcr above 
Tc for shghtly under doped HTSC. T^r is the T where 
p(T) deviates upward exponentially, which has been well 
explained (l|,|2,(3l via Ej in Eq. (|^ . 

8. Conclusions 

In conclusion, iFDS based electrical resistivity (with e- 
e scattering rate only) and heat conductivity (with both 
e-e and e-ph scattering rate) models have been utilized 
to tackle the T dependence of Lorentz ratio in optimally 
doped YBa2Cu307. The computed C{T) with experi- 
mental K overwhelmingly suggests that Bloch-Griineisen 
formula or the inclusion of e-ph scattering in the electrical 
resistivity is not suitable, at least for YBa2Cu307. On 



FIG. 1; Experimental 7?^''' (T) and pab{T) (inset) data points 
for YBa2Cu307_i single crystal (Al) have been fitted using 
Eqs. (I30II and 12911 respectively. The former equation is com- 
puted with two m*s namely, 50mo and 73mo while the resis- 
tivity is calculated with A = 7.3 x 10~^ 57-cm. 



FIG. 2: a) Shows the BG resistivity, p{T) plots above 90 
K for Qd = 350, 300 and 200 K. Whereas, b) depicts the 
theoretical plots for the BG Lorentz ratio, £bg above 90 K 
with experimental heat conductivity (k) using Eq. ^ with 
the Debye T, Qd = 350, 300 and 200 K. The calculated £(T) 
with Eq. 13111 using experimental k is also plotted with o 
in b). The theoretical solid line in b) satisfies iFDS based 
models namely, Eqs. (HHJ, (HU and ^ with a = 3. The 
experimental plots indicated with • is obtained from the data 
combined from Leridon et al. and Sutherland et al. G^. 

the other hand, e-ph scattering contributes significantly 
in heat conductivity that eventually gives a reasonably 
acceptable picture for the experimental heat conductivity 
and Lorentz ratio. Additionally, the spin Pseudogap phe- 
nomenon have been omitted throughout so as to avoid its 
inconclusive interpretations. Apart from that, the mag- 
nitudes of the T-independent scattering rate constant, ef- 
fective mass and the charge carriers density are all in the 
acceptable range, complying with other optimally doped 
YBCO single crystals as computed previously. 
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